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Using a sample of 471 million BE events collected with the BABAR detector, we study the sum of 
seven exclusive final states B -^ ^s(d)75 where Xs(^) is a strange (non-strange) hadronic system with 
a mass of up to 2.0 GeV/c . After correcting for unobserved decay modes, we obtain a branching 
fraction for h ^r d'y of (9.2±2.0(stai.)±2.3(sj/st.)) x 10~^ in this mass range, and a branching fraction 
for fe ^> S7 of (23.0 ± 0.8(stai.) ± 3.0(syst.)) x 10"^ in the same mass range. We find %ff^'ll^\ = 
0.040 ± 0.009(siaf.) ± 0.010(sysi.), from which we determine \Vtd/Vts\ = 0.199 ± 0.022(stoi.) ± 
0.024(syst.)±0.002(t/i.). 

PACS numbers: 

The decays b ^ dj and b ^ sj are flavor-changing Standard Model (SM). The leading-order processes are 
neutral current processes forbidden at tree level in the one-loop electroweak penguin diagrams, for which the 

top quark is the dominant virtual particle. In theories 
beyond the SM, new virtual particles may appear in the 

loop, which could lead to measurable effects on experi- 

*Now at Temple University, Philadelphia, Pennsylvania 19122, mental observables SUch aS branching fractions and CP 

y°^ asymmetries [l] . In the SM the inclusive rate for b ^- d'y 

T Also with Universita di Perugia, Dipartimento di Fisica, Perugia, • j i,. .?,, i r , \ir /Tr \2 



is suppressed relative to & — >■ 57 by a factor \Vtd/Vt 



ts\ 



Italy 

tAlso with Universita di Roma La Sapienza, 1-00185 Roma, Italy where Vtd and Vts are Cabibbo-Kobayashi-Maskawa ma- 

§Now at University of South Alabama, Mobile, Alabama 36688, trix elements. Measurements of | Vtd/ VtJ using the exclu- 

USA sive modes B — >• (p, uj)j and B -^ K*j [3,[1| are now well- 

^Also with Universita di Sassari, Sassari, Italy established, with theoretical uncertainties of 7% from 



TABLE I: The reconstructed decay modes. Charge conjugate 
states are imphed throughout this paper. 



B^Xaj 




B^Xsj 


B" -^ n+Tv-'j 




B" -^ K+n-'j 


B+ -^ TT + TT^ 




B+ -^ K+TT^-y 


B+ ^7r+7r"7r+7 




B+ ^ K+TT+Tv-y 


B" -^ TT + TV-TT^ 




B° -^ K+TY~TV°J 


B° -> TT + TV-TT + TV- 


'1 


B° ^ K+n-Tv+n-y 


B+ -,Tr + 7V-7T + 7v' 


'7 


B+ -^ A'+7r-7r+7r°7 


B+ -^ TT + TJ-y 




B+ -^ K+rjy 



weak annihilation and hadronic form factors [2]. This 
ratio can also be obtained from the Bd and Bs mixing 
frequencies [3| . It is important to confirm the consistency 
of the two methods of determining |Vtd/Vts|, since new 
physics effects would enter in different ways in mixing 
and radiative decays. A measurement of the branching 
fractions of inclusive b ^ d'j relative to 5 ^- 37 would de- 
termine |Vtd/Vts| with reduced theoretical uncertainties 
compared to that from exclusive modes [6j. 

This letter supersedes of Q, and presents the first 
significant observation of the 6 — ^ ^7 transition in the 
hadronic mass range AI{Xd) > l.OGeV/c^, resulting 
in a significant improvement in the determination of 
|Vtd/Vts| via the ratio of inclusive widths. Inclusive 
b -^ s"f and b -^ dj rates are extrapolated from the 
measurements of the partial decay rates to seven ex- 
clusive final states (see Table H]) in the hadronic mass 
ranges 0.5 < M{Xd) < l.OGeV/c^ (low mass, contain- 
ing the previously measured K*, p and lo resonances) 
and 1.0 < M{Xd) < 2.0GeV/c2 (high mass). We com- 
bine these measurements and make a model-dependent 
extrapolation to higher hadronic mass to obtain an in- 
clusive branching fraction [B) for b — >■ (5,^)7. These 
measurements use the full dataset of 471 x 10^ BB pairs 
collected at the T{AS) resonance at the PEP-II B factory 
with the BaBAR detector 0. 

High energy photons are reconstructed from an iso- 
lated energy cluster in the barrel of the calorimeter, with 
shape consistent with a single photon, and energy 1.15 < 
E* < 3.50 GeV, where * denotes the center-of-mass (CM) 
frame. We remove photons that can form a tt" (77) candi- 
date in association with another photon of energy greater 
than 30 (250) MeV if the two-photon invariant mass is in 
the range 110 < m^j < 160 (520 < nij^ < 560) MeV/c^ 
for the low mass region and 95 < m^y < 155 (530 < 
m-y-y < 565) MeV/c^ for the high mass region. 

Charged pion and kaon candidates are selected from 
well-reconstructed tracks. We use a pion selection algo- 
rithm to differentiate pions from kaons, with a typical 
selection efficiency of 95% and kaon mis-identification 
rate of 4%. Kaons are identified as tracks failing the 
pion selection criteria. We reconstruct n^{r]) candidates 
from pairs of photons of minimum energy 20 MeV with 
an invariant mass 115 < rrijj < 150 (470 < m-yy < 
620) MeV/c^ . We require all pion, 77 and kaon candidates 



to have a momentum in the laboratory frame greater than 
600 (425) MeV/c in the low (high) mass region. 

The selected pion, 77, kaon and high-energy photon can- 
didates are combined to form B meson candidates con- 
sistent with one of the seven decay modes. The charged 
particles are combined to form a common vertex with a 
X^ probability greater than 1%. We use the kinematic 
variables AE ^ Eg - E^^^^, where Eg is the energy of 
the B meson candidate and -Eboam i^ ^^^ beam energy, 
and rriES = V-^bcam ~ Pb^ ' where pg is the momentum 
of the B candidate. We consider candidates in the range 
-0.3 <AE< 0.2 GeV and tties > 5.22GeV/c2. 

Contributions from continuum processes (e+e^ — >■ qq, 
with q — u,d,s,c) are reduced by considering only events 
for which the ratio R2 of second-to-zeroth order Fox- 
Wolfram moments is less than 0.98. To further dis- 
criminate between the jet-like continuum background 
and the more spherically symmetric signal events, we 
compute the angle 9^ between the photon momentum 
and the thrust axis of the rest of the event (ROE) and re- 
quire I cos(6'5-)| < 0.8. The ROE is defined as all charged 
tracks and neutral energy deposits that are not used to 
reconstruct the B candidate. 

Ten other event shape variables that distinguish be- 
tween signal and continuum events are combined in a 
neural network (NN). These include the ratio i?2i which 
is i?2 is calculated in the frame recoiling against the pho- 
ton momentum, the B meson production angle with re- 
spect to the beam axis in the CM frame, 8*g, and the 
L-moments [T^] of the ROE with respect to either the 
thrust axis of the ROE or the direction of the high en- 
ergy photon. Differences in Icpton, pion and kaon pro- 
duction between background and B decays are exploited 
by including several flavor-tagging variables applied to 
the ROE fllj. Using the NN output, we reject 99% of 
continuum background while preserving 25% of signal de- 
cays 

After all selections are applied, there remain events 
with more than one B candidate. In these events the 
candidate with the reconstructed tt" or ry mass closest to 
nominal is retained. Where there is no tt*^ or 7/ we retain 
the candidate with the highest vertex y^ probability. 

The signal yields in the data for the sum of the seven 
decay modes are determined from two-dimensional ex- 
tended maximum likelihood fits to the AE and tt^es 
distributions. We consider the following contributions: 
signal, combinatorial backgrounds from continuum pro- 
cesses, backgrounds from other B decays, and cross-feed 
from mis-reconstructed B -^ X^ decays. The fits to 
B -^ Xd"f events contain components from misidentified 
b ^ s"f decays, and we neglect the small b ^ dj back- 
ground in the fits to B —>■ Xg'y events. 

Each contribution is modeled by a probability density 
function (PDF) that is determined from Monte Carlo 
(MC) simulated events unless otherwise specified. For 
the misidentified signal cross-feed components, we use a 
binned two-dimensional PDF to account for correlations. 
All the other PDFs are products of one-dimensional func- 



tions of Ai? and tties- For signal, the mp. s spectrum is 
described by a Crystal Ball function [lj|, and AE by 
a Cruijff function [1J|. The parameters of these func- 
tions are determined from the fit to the high-statistics 
B — >■ Xsj data sample. We use these fitted values to fix 
the signal shape in the fits to B ^ Xd'y events. 

The remaining B backgrounds contain a small compo- 
nent that peaks in tties but not AE, which is modeled by 
a Gaussian distribution in tties- Continuum and other 
non-peaking backgrounds are described by an ARGUS 
shape 14] in tties ^nd a second-order polynomial in AE. 

We perform separate fits for B — >• X^j and B — j- Xs'y 
in each of the hadronic mass ranges 0.5-1.0 GeV/c^ and 
1.0-2.0 GeV/c^. For each of the four fits, we combine 
the component PDFs and fit for the signal, generic B 
and continuum yields, the ARGUS and two polynomial 
shape parameters. We scale the cross-feed contributions 
proportionally to the fitted signal yield, re-fit and iterate 
until the procedure converges. Projections of toes and 
AE from fits to data for B -^ Xs'f and B -^ Xd'y are 
shown in the high mass regions in Figure [T] Table HIl 
gives the signal yields, efficiencies (after corrections for 
systematic effects) and partial branching fractions (VB) . 
We calculate VB using VB{B -^ X7) = Ns/{2 eNgs), 
where N^q is the number of BB pairs in the data sample. 
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FIG. 1: Projections of AE with 5.275 < rriES < 5.286 GeV/c^ 
for (a) B — >■ Xsj and (b) B — >■ Xd^, and of tties with 
-0.1 < AE < 0.05 GeV for (c) B ^ X.j and (d) B -^ Xa'y 
in the mass range 1.0-2.0 GeV/ c^. Data points are compared 
with the sum of all the fit contributions (solid line). The 
jagged line is an artifact of the fit projection over the sum of 
several binned histograms. The dashed line shows the signal 
component. 



We have investigated a number of sources of systematic 
uncertainty in the measurement of the partial branching 



fractions, some of which are common to both B — >■ Xdj 
and B -^ X^j and cancel in the ratio of branching frac- 
tions (see Table Hill those that do not cancel in the ratio 
are marked by an asterisk). Uncertainties in tracking ef- 
ficiency, particle identification, 7 and tt" reconstruction, 
and the n'^/r] veto have been evaluated using indepen- 
dent control samples of data and MC simulated events, 
and incorporated into our analysis. Uncertainty due to 
the NN selection has been evaluated by comparing the ef- 
ficiency of the selection in data and MC for the B — ^ Xgj 
events, which are relatively free of background, assuming 
that potential discrepancies between data and MC are 
the same for the B —^ Xdj sample. The means and 
widths of the signal PDF are varied within the range 
allowed by the fit to the B — > Xgj data, accounting 
for correlations. Other PDF parameters are also varied 
within the la limits determined from the fit to MC. We 
vary the 5 — )■ 57 background in the fit to i? — >■ Xdj by 
the statistical uncertainty on our measurement of those 
decays. The signal cross-feed originating from our mea- 
sured channels is varied by the statistical uncertainty on 
our measurement; other signal cross- feed backgrounds by 
±50%. An additional uncertainty on the efficiency arises 
from the fragmentation of the hadronic system among 
the measured final states. For B -^ Xgj the uncertainty 
is constrained by the errors on the measured data; for 
B — >■ Xdj an estimate is obtained from the difference be- 
tween the default phase-space fragmentation (see below) 
and a re-weighting using the measured data/MC differ- 
ences in S — > Xs7. 

To obtain inclusive B{b -^ 57) and B{b -^ ^7) we need 
to correct the partial B values in Table |ll] for the frac- 
tions of missing final states. After correcting for the 50% 
of missing decay modes with neutral kaons, the low mass 
B -^ Xsj measurement is found to be consistent with 
previous measurements of the rate for B — >■ K*j [l^. For 
the low mass B -^ Xdj region, we correct for the small 
amount of non-reconstructed uj final states (for example, 
u! — > 7r*'7), and find a partial branching fraction consis- 
tent with previous measurements oi B{B — >■ {p,uj)^) [15| . 
We assume that non-resonant decays do not contribute 
in this region. 

In the high mass region, the missing fractions depend 
on the fragmentation of the hadronic system and are ex- 
pected to be different for Xd and Xg- In our signal MC, 
fragmentation is modeled by selecting an array of final- 
state particles and resonances according to the phase- 
space probability of the final state, as implemented by 
JETSET [H]. We further constrain the distribution of 
Xg final states to that observed for our seven decay modes 
as well as the distributions of a number of other states 
measured in [l3|. According to this MC model we re- 
construct 43% of the total width in 6 — >■ ^7 , and 36% 
in 6 — > S7 . A further 37% of the width of 6 — > 37 is 
constrained by the isospin relation between charged and 
neutral kaon decays. We explore the uncertainty in the 
correction for missing modes by considering several al- 
ternative models: replacing 50% of 6 ^- 57 and b ^- dj 



TABLE II: Signal yields (Ns), efficiencies (e), partial branching fractions {VB), inclusive branching fractions {B) and the ratio 
of inclusive branching fractions for the measured decay modes. The first error is statistical and second is systematic (including 
an error from extrapolation to missing decay modes, for the inclusive B). 

M(X,)0.5-1.0 M(Xd)0.5 - 1.0 Af (X^)l.O - 2.0 M(Xd)1.0 - 2.0 M(X,)0.5-2.0 M(Xd)0.5 - 2.0 



Ns 

e 

VB{xW-'^) 

S(xlO-^) 

Bjb^dj) 

B(h->S7) 



804 ± 33 

4.5% 
19±1±1 
38 ± 2 ± 2 



35 ±9 
3.1% 

1.2 ±0.3 ±0.1 

1.3 ±0.3 ±0.1 



990 ± 42 

1.6% 

66 ± 3 ± 6 

192 ± 8 ± 29 



56 ±14 

1.9% 

3.2 ±0.8 ±0.5 

7.9 ±2.0 ±2.2 



230±8±30 9.2 ±2.0 ±2.3 



0.033 ±0.009 ±0.003 



0.040 ±0.009 ±0.010 



TABLE III: Systematic errors on the measured partial and 
inclusive branching fractions B. Systematic errors that do 
not cancel in the ratio of rates are marked with (*). 



Systematic 


M{Xs) 


M{Xa) 


Error Source 


0.5-1.0 


1.0-2.0 


0.5-1.0 


1.0-2.0 


Track selection 


0.3% 


0.4% 


0.3% 


0.4% 


Photon reconstruction 


1.8% 


1.8% 


1.8% 


1.8% 


tt"/?? reconstruction 


0.9% 


1.1% 


1.4% 


1.6% 


Neural network 


1.1% 


4.9% 


1.1% 


4.9% 


B counting 


0.6% 


0.6% 


0.6% 


0.6% 


PID (*) 


2.0% 


2.0% 


2.0% 


2.0% 


Fit bias (*) 


0.1% 


0.9% 


4.9% 


6.5% 


PDF shapes (*) 


2.3% 


0.6% 


3.7% 


3.4% 


Histogram binning (*) 


0.8% 


0.2% 


1.8% 


1.8% 


Background (*) 


0.8% 


1.2% 


5.9% 


7.0% 


Fragmentation (*) 


- 


3.3% 


- 


5.1% 


Signal model 


- 


5.8% 


- 


6.0% 


Error on partial B 


4.0% 


9.0% 


9.3% 


14.2% 


Missing > 5 body 




9.6% 




18.2% 


Other missing states 




7.5% 




15.3% 


Spectrum Model 




1.8% 




1.6% 


Error on inclusive B 


4.0% 


15.2% 


9.3% 


27.7% 



hadronic final states with a mix of resonances; varying 
6 — !■ S7 fragmentation constraints within their statisti- 
cal uncertainties; and setting the 6 — >■ ^7 fragmentation 
rates to those of their corresponding & ^ 57 states. The 
resulting missing fractions vary by up to 50(40)% rela- 
tive to the nominal model in i? — ;> Xsj{B -^ Xdj). We 
therefore independently vary final states with > 5 stable 
hadrons, or with > 27r° or rj mesons, by ±50(40)%. 

Results for the corrected B values are shown in Ta- 
ble HH Adding the two mass regions, taking into account 
a partial cancellation of the missing fraction uncertain- 
ties in the ratio oi b —i' dj to b —^ sj , we find B{b -> 
d-i)/B{b -^ 57) = 0.040 ± 0.009(stai.) ± 0.010(s2/si.) in 
the mass range M{X) < 2.0GeV/c^ 

We correct for the unmeasured region M{X) > 
2.0 GeV/c^ using the spectral shape from Kagan- 
Neubert [18'| with the kinetic parameters {mb,fJ,'^) = 



(4.65 ± 0.05, -0.52 ± 0.08) extracted from fits of 6 -J> 57 
and b -^ civ data [10], yielding corrections of 1.66±0.03; 
the spectra for b ^f s^ and b ^- d'y are expected to be 
almost identical. 

Conversion of the ratio of inclusive branching fractions 
to the ratio iVt^/^sl is done according to [a], which re- 
quires the Wolfenstein parameters p and fj as input. How- 
ever, since the world average of these quantities relies on 
previous measurements of \Vtd/Vts \ we instead re-express 
p and ?7 in terms of the world average of the indepen- 
dent CKM angle /3 15]. This procedure yields a value of 
\Vtd/Vts\ = 0.199±0.022(stai.)±0.024(sysi.)±0.002(i/i.), 
compatible and competitive with more mo del- dependent 
determinations from the measurement of the exclusive 
modes B -^ (p, 1:^)7 and B -^ K*j [1, [!]■ 

In summary, we have measured the inclusive 6 — )► 57 
and b —>■ d^ transition rates using a sum of seven fi- 
nal states in the hadronic mass range up to 2.0GeV/c^, 
making the first significant observation of the b ^ d'j 
transition in the region above l.OGeV/c^. The value of 
IVtd/VtsI derived from these measurements has an exper- 
imental uncertainty approaching that from the measure- 
ment of exclusive decays B -^ {p, ll!)j and B -^ K*j, but 
a significantly smaller theoretical uncertainty. 
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